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Abstract 

The understanding of the detailed procedures of the geodetic technique for structural deformation 
monitoring and analysis is required to enable a reliable and cost-effective monitoring system for the 
monitoring of engineering structures, be developed. To do this, the different stages that are involved in 
deformation studies using the geodetic method are explained in detail. As a result, this paper presents 
detailed geodetic technique procedures for structural deformation monitoring and analysis. Data 
collection, data processing, result presentation and analysis stages are discussed in detail. 
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1 Introduction 

The incessant collapse of engineering structures such as bridges, viaducts, dams, high rise buildings and 
tunnels in various parts of the world has resulted in their (the structures) monitoring for safety purpose. 
Structural, as well as deformation monitoring, is done to detect movement that could lead to collapse and 
to allow for sufficient warning to successfully evacuate the structure. As it is well known, engineering 
structures are subject of deformation due to factors such as changes of groundwater level, changes in the 
bedrock, increase or decrease of weight, changes of the material properties as a result of changes in 
temperature, ageing or outside influences, tectonic phenomena and human activities. 

The importance of deformation monitoring cannot be underestimated. Several cases of collapsed 
structures in Nigeria and various parts of the world have been reported. In 2017, a bungalow with a shop 
extension collapsed in Lagos, in 2018 a seven-story building collapsed in Port-Harcourt. Also in 2018, a 
building marked for demolition in Miami, United States, collapsed. In Italy, Morandi Bridge in Genoa, 
collapsed, which led to the death of 43 people. Abdullahi and Yelwa (2016) also gave several instances 
of collapsed structures as those of the Ojiani dam in Akoko-Edo, in Edo State, the Tiga dam in Kano 
State some years ago which led to the loss of properties and displacements of thousands of rural 
inhabitants and the collapse of buildings in Abuja, Lagos and Ibadan. Others are Highland Towers 
collapse in Malaysia, Savar building collapse in Bangladesh, Omsk building collapse in Russia, Royal 
Plaza Hotel collapse in Thailand, Grace Divine School collapse in Haiti, etc. All these could have been 
averted if there had been a program for monitoring the stability or deformation of these structures, 
according to Abdullahi and Yelwa (2016). This was as a result of undermining the importance of 
deformation survey, and not knowing and understanding the procedures for the development of a reliable 
and cost-effective monitoring system for the monitoring of engineering structures. 

Most engineering structures are monitored using the geodetic method such as the digital level for vertical 
displacement and the total station and Global Positioning System (GPS) for horizontal displacement 
detection. The selection of the method of measurements depends upon the accuracy requirements for the 
survey (Abdel-Gawad et al, 2014). The implication of using deformation monitoring method that is not 
suitable in term of accuracy for the monitoring of a particular engineering structure has a significant 
effect on the magnitudes and directions of the determined deformations (movements) (Eteje, 2019). 
According to Engineer Manual (2018), surveying accuracy specifications are meant to ensure detection 
of a given amount of movement under normal operating condition. 

The analysis of the deformations, as well as movements of the structures, is crucial to decide whether a 
significant movement has taken place between observations epochs. This is because the said 


7 | P a g e 
www.iiste.org 



International Journal of Scientific and Technological Research 
ISSN 2422-8702 (Online), DOI: 10.7176/JSTR/6-07-02 
Vol.6, No.7, 2020 


www.iiste.org 


W II 

iisTe 


deformations are often of the same magnitudes as the observations errors. To do this, some statistical 
analyses have to be carried out to confirm if a significant movement has taken place or not. 

To develop a reliable and cost-effective monitoring system for the monitoring of engineering structures, 
the following procedures are considered: the deformation monitoring scheme must consist of discrete 
object monitoring points located on the surface of the structure where the movement of the entire 
structure will be determined, a set of stable control points established which defines the reference frame 
for the movements (Ebeling, 2014), geodetic observations collected at discrete time intervals, or epochs 
that describe the relative geometry between the monitoring and control/reference points, a set of 
coordinates and their accuracy estimated from these observations that describe the state of the monitoring 
network at each epoch. If data from multiple epochs are available, deformations that occurred between 
these epochs can be derived from two epoch analysis (Ebeling, 2014). 

To understand the detailed procedures for the development of a reliable and cost-effective monitoring 
system for the monitoring and analysis of engineering structures for safety purpose, this paper presents 
detailed geodetic technique procedures for structural deformation monitoring and analysis. 

2 The Detailed Procedures for Structural Deformation Monitoring and Analysis 

2. 1 Monitoring Points Located on the Surface of the Structure 

The first step is to select/choose suitable points on the structure where studs would be placed. These 
points must be selected such that they can be observed from not less than two reference/control points. 
This is necessary because of the least squares adjustment of the observations. The studs must be placed 
fixed such that any of their movement would be the displacement of the structure itself. 

2.2 Establishment of Control/Reference Points 

A set of control/reference points has to be established on stable grounds or platforms to define the 
reference frame for the movements. The control points are usually monumented and assumed to be stable 
throughout the monitoring period. They are also selected where they will not be disturbed and where 
observations can easily be taken from them to not less than two monitoring points. The observation of 
the control stations should be done accurately. It should be carried out relative to a set of first-order 
controls using either GNSS observation or a high precision total station. 

2.3 Geodetic Observations Collection at Discrete-Time Intervals or Epochs 

Having selected the monitoring points on the structure and established the reference stations around the 
structure, the next step is to carry out observations from the reference stations to the monitoring points 
at a specified interval. The interval may either be one month, three months, six months or one year 
depending on the purpose of the survey and the magnitude and direction of the expected 
deformation/movement. In each epoch observations, each of the monitoring points must be observed 
from more than one reference station to allow for redundant observations, as well as least squares 
adjustment of the epoch observations (see Figures la and lb). 



Fig. la: Taking Total Station Observations from the Reference Stations to the Monitoring Points 
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Fig. lb: GNSS Observation of the Monitoring Points Relative to the Control Stations 
2.4 Processing of Each Epoch Observations 

Each epoch observations are processed and adjusted to obtain: the adjusted coordinates and heights of 
the monitoring points, the standard errors of the positions and heights of the monitoring points, the error 
ellipses of each of the processed positions, heights confidence intervals, the a posteriori variance and a 
posteriori standard error of the observations, the trace of the variance-covariance matrix, Test of the 
Variance of a Probability Distribution (Chi-Square Test), Redundancy Number, Outlier Detection, 
Residual Distribution (Goodness of Fit Test) and Tau Statistics. The adjustment is carried out using the 
observation equation method of the least squares technique. 

2.4.1 Observation Equation Method of Least Squares Adjustment 

The observation equation method of least squares adjustment technique is applied in deformation studies 
to adjust each epoch observations to enable the precision and accuracy of the adjusted positions and 
heights, as well as those of the observations, be obtained. In the observation equation method, there is a 
functional relationship between every measurement and certain parameters. Every measurement 
(observation) can be related to a certain number of parameters by an equation and the number of equations 
that can be written, therefore, is exactly equal to the number of measurements used in the adjustment. A 
solution to these equations is obtained by minimizing the sum of the (weighted) squares of residuals 
(least squares principle) (equation 1). This Teast squares’ solution results in direct determination of 
parameters. 


n 


</) = v\ + + v 2 + ...+v 2 = £v 2 - Minimum 


( 1 ) 


Ono et al. (2014) gave the following functional relationship between the adjusted observations and the 
adjusted parameters as: 


( 2 ) 


La = F(X a ) 


Where, L a = adjusted vector of observations and X a = adjusted station coordinates. Equation (2) is a linear 
function and the general observation equation model was obtained. 

To make the matrix expression for performing least squares adjustment, an analogy will be made with 
the systematic procedures. The system of observation equations is presented by matrix notation as 
(Mishima and Endo 2002, Ono et al, 2018, Eteje and Oduyebo, 2018) 


V=AX-L 


( 3 ) 


Where, 
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Estimated parameter 

X=0v)A) (5) 

Where, 

N = ( A T WA ) = Normal Matrix 

t = (a t wl) 

AT 1 = (a t WaY = Q „ 

X = (a' WA) ] (a' WL) = Estimate 
W = Weighted Matrix 

2.4.2 Observation Equation 

The application of least squares adjustment technique for the adjustment of observations in deformation 
studies requires that observation equations should be formed or written. The number of observation 
equations must be equal to the number of observations collected. Here, observation equations for GPS 
baseline vectors for GNSS method, distance and azimuth observation equations for total station method 
and levelling circuit observation equations are detailed. 

2.4.2.1 Least Squares Adjustment Model for GPS Baseline Vectors 

Ghilani and Wolf (2008) explained that from the known X, Y and Z coordinates of stations A and B , and 
the observed AX, AY and AZ components, coordinates of new stations C, D, E and F can be calculated. 
An adjustment is necessary because redundant observation exists. In applying least squares to the 
problem, observation equations are written that relate the unknown adjusted coordinates of the new 
stations C, D, E, and F to the observed AX, AY and AZ values and their residual errors. The following 
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observation equations as given by Ghilani and Wolf (2008) are written for the first two baselines (AC 
and DC) as 
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Ghilani and Wolf (2008) expressed the observations in equation (6) in matrix form as 
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2. 4.2.2 Distance Observation Equation 

According to Githumbi (2014), the distance observation equation is 

s ij =s: h + eo 

The distance S 12 between two points i and j can be computed as 

s»°=|x,-x,Mi 'i-rjf 1 


(7) 


( 8 ) 


Equation (8) is not linear with respect to the unknown parameters (X’s and Y’s) and therefore cannot be 
subjected to a least squares process before it is linearized. The linearization is done by the Taylor series 
expansion of functions only up to the linear terms. Thus, we now get derivatives of .S’, ; with respect to X-, 
, Y : , Xj, Y j; as 
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Where S ° is the distance computed from approximate coordinates, X °, Y °, X °, X ° . 

y 1 1 j j 


2.4.2.3 Azimuth Observation Equation 

According to Mikhail and Gracie (1981), the azimuth of a line from point i to point j (in that specific 
direction) is defined as the horizontal angle measured clockwise from the reference meridian to the line. 
The direction Ay between two points i and j can be computed as (Githumbi, 2014): 


A u = tan- 


x j~ x t 

Y j~ Y i 


( 10 ) 


Mikhail and Gracie (1981) explained that the azimuth can have a value anywhere between 0° and 360°. 
It is important to make sure that the quadrant of Ay is correctly determined. This means that the sign of 

both the numerator X - — X [ and the denominator Y . — F of the tangent function must be taken into 
account. They further stated that if a computational aid, such as hand calculator or computer, is used and 
it yields only a positive or negative acute angle a , an appropriate constant C must be added to a to 
obtain the azimuth Ay. That is, 

Aij =a ij+ C (11) 


Ghilani and Wolf (2008) also stated that if the azimuth of the line is between 0° and 90°, the value of C 
is 0°. If the azimuth of the line is between 90° and 270°, C is 180° and if the line's azimuth is between 
270° and 360°, C is 360°. The linearized form of equation (11) is given in Githumbi (2014) as: 
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2.4.2.4 Least Squares Adjustment of Levelling Circuits 

According to Ghilani and Wolf (2008), in adjusting level networks, the observed difference in elevation 
for each course is treated as one observation containing a single random error. Observation equations 
are written that relate these observed elevation differences and their residual errors to the unknown 
elevations of the benchmarks involved. These can then be processed using the least squares adjustment 
technique given in section 2.4.1 to obtain the adjusted values for the benchmarks and their standard 
deviations. Ghilani and Wolf (2008) posited that weights in differential levelling are inversely 
proportional to course lengths. The course length must be in kilometres. Figure 2 shows a levelling 
network for the least squares adjustment. 



Fig. 2: Levelling Network for Least Squares Adjustment. 
Source: Ghilani and Wolf (2008) 
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2.4.3 Weight of Observation 

Weight is the worth or reliability of one measurement relative to a standard or another measurement. 
Measurements with higher precision (smaller standard deviation) should be assigned higher weights and 
measurements with higher weights should receive smaller corrections after an adjustment. They are 
relative, and therefore, determined by comparing with another measurement. If the observations have 
been made with the same precision in each case, it would seem fairly obvious to allot weights 
proportional to the number of observations in each set. Allan, (2013) explained that weights are assigned 
to observations so that the weight of observation is proportional to the inverse expected (prior) variance 
of that observation, 

W = 2- (13) 

a 

When a reference variance or variance factor ( <J Q ) is used, weight is given as: 

CTo 

W = — (14) 

<7 


The weight of a mean value computed from repeated measurements is proportional to the number of 
repetitions. The weighted matrix of uncorrelated measurements is a diagonal matrix given as: 
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(15) 


2. 4.4 Degree of Freedom / Redundant Observation 

Redundant observations are additional measurements which are not required for computing the desired 
unknown quantities, the parameters. Network redundancy, also known as the number of degree of 
freedom, is defined as the number of observations minus the number of unknown parameters (Equation 
16). The redundant of the observations is useful tools to determine the unknown parameters using the 
least squares adjustment techniques in geodetic networks (Mohammad, 2016). When there are redundant 
observations available, traditional coordinate geometry formulae can yield many possible coordinate 
solutions, depending on which observations are used as input. This inconsistency is remedied by least 
squares adjustment (Bird, 2009). 

r = m-n (16) 

Where, 

r = Degree of freedom or redundant observations 
m = Number of observations 
n = Number of unknown 


2.4.5 Measure of Accuracy of Adjusted Parameters 

The accuracy information about the estimated parameters is all contained in the variance-covariance 
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The accuracy information of the adjusted parameters that contain in the variance-covariance matrix 
include the: 

1. Standard errors of the adjusted positions and heights of the monitoring points. 

2. Error ellipses of the adjusted positions. 


13 | P a g e 

www.iiste.org 







International Journal of Scientific and Technological Research 
ISSN 2422-8702 (Online), DOI: 10.7176/JSTR/6-07-02 
Vol.6, No.7, 2020 


www.iiste.org 


W II 

iisTe 


3. Heights confidence intervals. 

2.4.5.1 Standard Errors of the Adjusted Positions and Heights of the Monitoring Points 
The standard errors of the adjusted coordinates and heights of the monitoring points are the measure of 
the accuracy of the adjusted parameters. Eteje and Oduyebo (2018) gave the model for the computation 
of the standard errors of the adjusted parameters as: 



Where, 

Q nn is a diagonal element of the inverse of the normal matrix ( N ~ l ). 


2.4.5.2 Error Ellipses of the Adjusted Positions 

In surveying and geodesy terms, error ellipse simply is a graphical representation of the marginal standard 
deviation of the parameters and the covariance between them. Error ellipses are usually computed during 
adjustment of horizontal or three - dimensional network. They allow a convenient way for the 
interpretation of the directional station position accuracy. Error ellipses are also used in optimising a 
network. They can be computed for each adjusted station and between connected stations. The stations' 
error ellipses are computed from the a posteriori covariance matrix for each station. To compute the 
relative error ellipse between stations, the a posteriori covariance matrix for the AT station, the TO station 
and the a posteriori covariance matrix which correlates the two stations are used. Because the error ellipse 
is a two-dimensional statistical result, the 2D expansion factor is used to scale the semi-major and semi¬ 
minor axes to the chosen confidence level (in this case, 95%). 
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Mikhail and Gracie, (1981) further pointed out that the variances G , and G , are the eigenvalues of 


the covariance matrix of the random vector 


X 

Y 


Mikhail and Gracie, (1981) also stated that the probability of the position given by X and Y lies on or 
within the error ellipse is 
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Where, c = Expansion Factor (for 95%, c = 2.448). 

The probability P^CJ < C 2 ] is represented by the volume under the bivariate normal density surface 
within the region defined by the error ellipse. The semi-major axis and semi-minor axis confidence ellipse 
is determined by 
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Semi-major axis confidence ellipse = G 7 C (22) 

Semi minor axis confidence ellipse = G 7 C (23) 

2. 4.5.3 Height Confidence Interval 

The Height Confidence Interval is an estimate of the precision of the adjusted height (orthometric or 
ellipsoidal, depending on the network adjustment context), at each station and the relative precision of 
the height between connected stations. The station height confidence interval is computed from the a 
posteriori covariance matrix for each station. To compute the relative height confidence interval between 
connected stations, the a posteriori covariance matrix for the AT station, the TO station and the a 
posteriori covariance matrix which correlates the two stations are also used. Since the height confidence 
interval is a one-dimensional result, the ID expansion factor (1.959) is used to scale the height confidence 
interval to the chosen confidence level (95%). The height confidence interval is determined by 
multiplying the adjusted height standard deviation by the confidence level (say 95%) expansion factor (c 
= 1.959 for ID) (Columbus, 2009). 

2.4.6 A Posteriori Variance and a Posteriori Standard Error 

The a posteriori variance which is also called the variance of unit weight and the a posteriori standard 
error which is also known as the standard error of unit weight are respectively measures of precision and 
accuracy of a given set of observations. In other words, they are respectively the average, as well as the 
mean precision and accuracy of a given set of observations. The models for the computation of the a 
posteriori variance and a posteriori standard error as given in Ameh (2013) are: 

1. A Posteriori Variance ( G o ) 



V T WV _ V T WV 
m — n r 


(24) 


2 . 


A Posteriori Standard Error ( A- ) 
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cr 
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V T WV 


m — n 


V T WV 


(25) 


Where V, W and m — n = r as respectively given in equations (3), (8) and (9). 


2.4.7 Variance- Covariance of the Adjusted Observations 

The model for the computation of the variance-covariance of the adjusted observations ^ ) is given 


as: 


y = A y a t = <j 0 2 an ~ 1 a t 

t-ULa X-UXa 


(26) 


2.4.8 Trace of the Variance-Covariance Matrix 

The trace of a square matrix A is defined as the sum of the diagonal elements of A, written as tr A or 
sometimes tr (A), that is, 


trA = 


(27) 


If A is a covariance matrix, then tr A is the sum of all variances and can be interpreted as a measure of 
the overall accuracy of the associated vector of random variates (Caspary, 1988). If two sets of adjusted 
parameters are given, the one with the smaller trace has better accuracy, hence, the smaller the trace the 
better the accuracy. 

2.4.9 Test of the Variance of a Probability Distribution (Chi-Square Test) 

Bird (2009) stated that after adjustment, certain statistical hypothesis tests are done to check the quality 
of the solution and to try to detect the presence of outliers. The confidence level normally used for 
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network adjustments is 95%. The corresponding significance level of 5% means that there is a 5% chance 
of making a type I error (Bird, 2009). A type I error is when a correct null hypothesis is rejected. 

Under the assumption that a population is normally distributed, we can test the null hypothesis H a that 

the population variance is <J 2 against the alternative that it is not <J 2 , using the sample variance S 2 as 
the test statistics (Mikhail and Grade, 1981). 


H 0 is accepted if S 2 the specific value of S 2 calculated from the sample lies between 


y cr o 

A- a,r 
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and 
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(28) 


otherwise, H a is rejected. 

Where, 

(J 2 = a priori variance 

S 2 = S' 2 = a posteriori variance 
r = degree of freedom or redundant observation 
1 — CC is the probability that H a is accepted when it is true. Thus, 1 — OC —b —a 

2.4.10 Redundancy Number 

The redundancy number for each adjusted observation is yet another statistic that may help to isolate 
problem areas in the network. For each adjusted observation in the network, a redundancy number is 
calculated. The redundancy number ranges from 0 to 1. Observations redundancy numbers are a measure 
of how close the variance of the residual is to the variance of the observation. If the redundancy number 
is close to 1, then the variance of the residual is close to the variance of the observation and the variance 
of the adjusted observation is close to zero. If the redundancy number is close to zero, then the variance 
of the residual is close to zero, and the variance of the adjusted observation is close to the variance of the 
observation (Leick, 1990). Intuitively it is expected that the variance of the residuals and the variance of 
the observations are close. When this is the case, the noise in the residuals equals that of the observations, 
and the adjusted observations are determined with high precision. The case where the redundancy number 
is close to 1 is preferred and it is said that the gain on the adjustment is high. If the redundancy number 
is close to zero, one expects the noise in the residuals to be small. 

The redundancy number of each adjusted observation ( ) and the average redundancy number of a group 

of adjusted observations ( r Av ) are respectively given as (Leick, 1990) 

r i = ViPi 

0 < r- <1 

n-R(A) 

r Av ~ 

n 

Where, 

q t = Diagonal element of the estimated residual cofactor matrix, Q v 
p. = Weight of the ith observation. 
n — R( A) = Degree of freedom or redundant observation. 

R(A) = Number of unknown parameters. 
fl = Number of observations. 

2.4.11 Outlier Detection 

According to Githumbi (2014), an outlier is defined as a residual, which according to some test rule is in 
contradiction to the reference. This operational definition of an outlier allows for a test strategy and a 


(29) 

(30) 
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clear statistical concept, but is only a relative definition depending on the selected risk level, the assumed 
distribution and the test procedure. Despite this fundamental difference between the definition of outliers 
and gross errors, it is naturally expected that the detected outliers are caused by gross errors. 

Whenever outlying residuals are detected, it is necessary to thoroughly check the records of the 
observations to find out if there are any gross errors. If no gross error can be found, the corresponding 
observation is detected and if necessary, re-measured. The approach to outlier detection is based on the 
Gauss Markov Model. The first step is usually a global test. If this global test fails, a procedure to flag 
erroneous measurements follows i.e. robust estimation. 

2.4.12 Residual Distribution (Goodness of Fit Test) 

The Standardized Residual for an observation is simply the ratio of the adjusted observation residual to 
the residual standard deviation that is, 

„ , , _ ., , Adjusted Observation Residual (v,) 

Standardized Residual (T t ) = —-- 

Residual Standard Deviation ( <j v ) 

The Standardized Residual test is one of the better indicators for removing possible outlier observations. 
If the Standardized Residual exceeds the Standardized Residual Rejection Constant (Tau Statistic), the 
Standardized Residual value is rejected. 

The Standardized Residual is a unitless quantity that allows both angular and linear observations to be 
compared directly. Probability theory suggests that random errors taken from large samples will tend to 
be normally distributed about their mean. Graphically, it can be represented by a bell-shaped curve, where 
the majority of the standardized residuals are expected to be clustered about the midpoint by combining 
a histogram and the numerical results from a Chi-squared Goodness of Fit” test. All the standardized 
residuals are grouped into classes about the theoretical mean (zero). Any standardized residual outside 
the leftmost or rightmost class is tabulated as out of range. Therefore, the closer the histogram 
approximates the bell-shaped curve, the closer the standardized residuals approach normality. 

The Chi-squared Goodness of Fit test is a quantitative measure as to the fit of the histogram to the normal 
distribution. Columbus (2009) stated that the graphical representation of the residual distribution should 
be considered in combination with the numerical Chi-square test results before concluding the residuals 
are or are not normally distributed. Columbus (2009) further stated that even if the test fails, the residuals 
may be considered normally distributed if they form a bell-shaped curve. 

2.4.13 Tau Statistics 

The Tau statistics are used to determine the critical value, as well as the residual rejection constant in 
residual distribution. Pope (1976) and Bird (2009) defined Tau statistics in terms of Student’s distribution 
as: 



(31) 

Where, 

T v a = Critical value from the Tau distribution at redundancy v and significance level Ot 

v = Degree of freedom 
OL = Significance level 

t v a = Critical value from the student's distribution at redundancy v and significance level CL . 
According to the rule of the test, H a is rejected if 7) is greater than T v a (T i > T va ) otherwise H a is 
accepted. Where T t is the standardized residual of each observation. 

2.5 Computation of Magnitude and Direction of Displacement/Movement 

Having processed each epoch observations to obtain the adjusted coordinates of the monitoring points 
and their respective standard errors, and the mean accuracy of the observations, the next step is to 
compute the magnitude and direction of the displacement/movement that takes place between 
observation epochs. The monitoring point displacement, d is computed as the difference in 
coordinates/heights between the measurement epochs as given by Eteje et al. (2018) is 
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Where, 

i i 

X c , y c = Coordinates of last epoch. 

P P 

X c , y c - Coordinates of the preceding epoch. 

And the horizontal displacement magnitude is computed as (Kaloop and Li, 2009) 

|d| = V^ 
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(32) 


(33) 


The vertical movement ( dH) is computed for each object (monitoring) point as (Eteje et al ., 2018): 


dH = V(4 -z*J = M 2 

Where, 

z\ - Height of last epoch. 

Z v c = Height of preceding epoch. 

The direction of movement (a) is computed as: 



(34) 


(35) 


2.6 Computation of 95% Confidence Ellipses/Intervals 

The 95% confidence ellipses/intervals are computed with the standard errors of the adjusted parameters 
and the 95% confidence level expansion factor, 1.96. In this, the error associated with the observed 
monitoring point in two different epochs is computed and multiplied by 1.96. The model for the 
computation of the 95% confidence ellipse is given by Beshr and Kaloop (2013) as 

e„ =1.96 faj +of (36) 


2 2 

Where, C y is the standard error squared in position for the (final) or most recent survey, cr. is the 

2 2 2 

standard error squared in position for the (initial) or reference survey. That is, CJj — <J^ and 

2 2 2 

cr. = <J iN + <J iE . This implies that, 

e n — 1 .96.^1 cr^ + cr^ + <J iN + o iE (37) 

Where, 

2 

<j^ = Standard error squared in northing of the final survey 

2 

<3j E = Standard error squared in easting of the final survey 

<J iN = Standard error squared in northing of the initial survey 

2 

<J iE = Standard error squared in easting of the initial survey 
2.7 Analysis 

This consists of the analysis of the adjusted observations which in turn comprises the analysis of the 
standard errors of the adjusted positions and heights, error ellipses of the adjusted coordinates, height 
confidence intervals, a posteriori standard error and a posteriori variance of the adjusted observations, 
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the trace of the variance-covariance matric, Chi-square test, redundancy number, outlier detection, the 
goodness of fit test and Tau statistics. It also consists of two epoch analysis. 

2.7 .1 Analysis of the Adjusted Observations 

The standard errors of the adjusted positions/coordinates and heights are analysed based on the computed 
values. It is a measure of the accuracy of the adjusted positions/heights. The smaller the computed values 
of the standard errors, the better the accuracy of the adjusted coordinates and heights. It is an indicator 
of the confidence we have in the adjusted parameters. If the computed values of the standard errors of 
the adjusted positions and heights are very small, then the confidence in the adjusted parameters is high 
and they are accepted. 

The error ellipses of the adjusted positions are used to present graphically the directional accuracy of the 
adjusted positions. It consists of the direction of the errors and the computed standard errors of the 
adjusted coordinates. When computed and plotted, it enables one to infer that the computed error is in a 
particular direction say north-east, north-west, south-east or south-west as the case may be (see Figure 
3). 



A- 


Fig. 3: Plot of Error Ellipses 


The height confidence interval is like the positions error ellipses, just that it is one dimension. It is also 
used to present graphically the accuracy of the adjusted heights. Since height is a vertical component, the 
heights confidence intervals are plotted vertically (see Figure 4). The smaller the plotted confidence 
intervals, the higher the accuracy of the adjusted heights. 


19 | P a g e 

www.iiste.org 




International Journal of Scientific and Technological Research 
ISSN 2422-8702 (Online), DOI: 10.7176/JSTR/6-07-02 
Vol.6, No.7, 2020 


www.iiste.org 


W II 

IIST| 



Fig. 4: Plot of Height Confidence Intervals 

The a posteriori standard error and a posteriori variance of the adjusted observations are respectively 
measures of the mean/average accuracy and precision of the adjusted epoch observations. They are 
computed with the residuals and the degree of freedom, as well as the number of redundant observations. 
If the values of the computed a posteriori standard error and a posteriori variance are very high, it implies 
that the accuracy and precision of the epoch observations are low. This means that the residuals, as well 
as the variations between the adjusted values and the observed values, are very large. 

The trace of the variance-covariance matric is another measure of accuracy. It is the sum of the elements 
of the principal diagonal of the variance-covariance matrix. The principal diagonal of the variance- 
covariance matrix contains the variations between the adjusted and the observed values. If the differences 
between the adjusted and the observed values are very small, it implies that the collected observations 
are very accurate. Their sum is also regarded as the overall accuracy of the adjusted observations. So, the 
smaller the trace of the variance-covariance matrix, the higher the accuracy of the collected observations. 
Chi-square test on variance factor is normally carried out to determine how the adjusted result matches 
the expected result. In this, a priori variance is used which depict the expected variation. In most cases, 
a priori variance of 1 is used. If the computed a posteriori variance falls within the range of the upper and 
lower limit as given in equation (28), the set of observations are accepted. But if it falls outside the limit, 
two things involve, it is either the a priori variance is too large or very small. If the a priori variance is 
too large, it implies that the accuracy of the collected observations is higher than the expected one. 
Therefore, the observations are accepted. But if on the other hand, the a priori variance for the survey is 
very small, it means that very high accuracy is expected. And when this accuracy is less than the expected 
one, the a posteriori variance will fall outside the range. 

Redundancy number is computed to determine the accuracy of individual observation. It ranges from 0 
to 1. If the redundancy number is close to 1, then the variance of the residual is close to the variance of 
the observation and the variance of the adjusted observation is close to zero. If the redundancy number 
is close to zero, then the variance of the residual is close to zero, and the variance of the adjusted 
observation is close to the variance of the observation (Leick, 1990). Here, if the redundancy number is 
close to 1 is better. In a group of observations where some of the redundancy numbers are close to 1 and 
others are close to 0 (zero), if the average redundancy number is close to 1, the set of observations are 
accepted. 

Outlier detection is done to determine the gross errors in a given set of observations. And this can be 
done using the goodness of fit test and Tau statistics. In the goodness of fit test, standardized residuals 
are plotted as a histogram, grouped into classes and combined with a normal distribution, as well as a 
bell-shaped curve. If all the histogram bars fall within the normal distribution curve, it means that the 
residuals are normally distributed and the observations that they represent are accepted. But if any of 
them falls outside the bell-shaped curve, the observations in that class are regarded as outliers as they 
contain gross errors and are rejected. This test is normally combined with the numerical Chi-square test 
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(see Figure 5). Columbus (2009) stated that even if the Chi-square test fails, the residuals may be 
considered normally distributed if they form a bell-shaped curve. Hence the observations are accepted. 



Fig. 5: Goodness of Fit Test 

In the Tau statistics, the computed standardized residuals are compared with the standardized residual 
rejection constant. If the computed standardized residual is less than the standardized residual rejection 
constant, the observation is accepted, otherwise rejected. Detected gross errors are removed from a given 
set of observations as their presence affects the overall accuracy of the observations. 

2.7.2 Deformation (Two Epoch) Analysis 

Deformation analysis is often regarded as two epoch analysis. It is done to determine if significant 
movement/displacement has taken place between any two observation epochs. To do this, the computed 
displacement for each monitoring point and its corresponding associated error as respectively given in 
sections 2.5 and 2.6 are compared. If the computed displacement is greater than the associated error, as 
well as the 95% confidence interval/ellipse, then a significant movement has taken place between the 
two measurement epochs. But if on the other hand, the computed displacement is less than the associated 
error, as well as the 95% confidence interval/ellipse then significant movement did not take place 
between the two measurement epochs. Therefore, the computed displacement is regarded as 
observation/measurement errors. 

3 Conclusion 

In conclusion, this paper has given the step by step procedures to consider developing a reliable and cost- 
effective geodetic monitoring system for the monitoring of engineering structures. The selection of 
monitoring points on the structure, establishment of reference stations around the structure, taking 
observations from the reference stations to the monitoring points at a specified interval, processing of 
each epoch observations, analysis of each epoch observations and two epochs analysis are discussed in 
detail. 
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